2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



1 



The cosmic-ray puzzle and the census of the interstellar medium: the 
Fermi LAT view of Cassiopeia, Cepheus and the Perseus arm. 

L. Tibaldo^'^ and I. A. Grenier^ 

on behalf of the Fermi LAT collaboration 

^ INFN - Sezione Di Padova & Dipartimento di Fisica "G. Galilei" - Universita di Padova, 
i -351 31 Padova, Italy 

^ Laboratoire AIM, CEA-IRFU/CNRS/Universite Paris Diderot, Service d'Astrophysique, CEA Saciay 
91191 Gifsur Yvette, France 



O 

(N 

00 
(N 

Oh 
O 

(N ■ 
> 

oo 

■ 

m ; 

On : 
O 



X 



Diffuse 7-ray emission arising from interactions of cosmic rays with the interstellar gas traces the densities of 
both of them throughout the Milky Way. We discuss the results obtained from the analysis of Fermi LAT 
observations in the region of Cassiopeia and Cepheus, towards the Perseus spiral arm. We find that the 7-ray 
emissivity of local gas is consistent with expectations based on the cosmic-ray spectra measured at the Earth. 
The emissivity decreases from the Gould Belt to the Perseus arm, but the measured gradient is flatter than 
predictions by a propagation model based on a cosmic-ray source distribution peaking in the inner Galaxy as 
suggested by pulsars. The Xco = ^(H2)/Wco conversion factor moderately increases by a factor ^ 2 from 
the Gould Belt to the Perseus arm. The presence of additional gas not properly traced by H I and CO surveys 
in the Gould Belt is suggested by the correlation between 7-rays and thermal emission from cold interstellar 
dust. 



1. INTRODUCTION 



More than a century after the discovery of cosmic 
rays (CRs), their origin is stih mysterious and their 
propagation in the interstehar space subject of long- 
standing debates. Interactions between CRs and in- 
terstehar gas give rise to 7 rays through 7r° decay and 
Bremsstrahlung radiation. In addition, CR electrons 
and positrons produce 7 rays through inverse Comp- 
ton scattering on the low energy interstellar radiation 
field. Thus, 7 rays can probe CR densities in the solar 
neighborhood, beyond the direct measurements in the 
solar system, as well as in remote parts of the Galaxy. 

On the other hand, the interstellar 7-ray emission 
usefully complements gas and dust tracers at other 
wavelengths to probe the interstellar medium (ISM) 
column densities. Two main radio/microwave tracers 
of the ISM have long been used in 7-ray astrophysics: 
the 21 cm line of atomic hydrogen. Hi, and the 2.6 
mm line of carbon monoxide, CO. There are, however, 
open issues in the derivation of gas column densities 
from these observations, and hints that their combi- 
nation does not provide an exhaustive census of the 
gas in the ISM. 

In this contribution, we add a few pieces to this 
intriguing puzzle with the use of Fermi Large Area 
Telescope (LAT) data for the second Galactic quad- 
rant, at 100° < / < 145°, -15° <b< 30°. This re- 
gion has been selected due to the presence of conspic- 
uous atomic and molecular complexes and the good 
kinematic separation between the different structures 
seen along the line of sight. Moving toward the outer 
Galaxy, the region encompasses complexes that are 
part of the Gould Belt (within ^ 300 pc from the So- 
lar System), of the local arm 1 kpc), of the Perseus 
arm (2.5—4 kpc) and of the outer spiral arm (> 5 kpc). 



The analysis is described in detail in [3| . We summa- 
rize below the main results and discuss their physical 
implications. 



2. GAMMA-RAY EMISSION FROM THE 
LOCAL INTERSTELLAR MEDIUM 

The 7-ray emissivity of atomic gas can be estimated 
by comparing 7-ray and radio data because the 21 cm 
line of the H i hyperfine transition allows us to derive 
hydrogen column densities, A^(Hi). The results pre- 
sented here are based on the often-used assumption (in 
7-ray astrophysics) of a uniform H i spin temperature 
of 125 K. There are, however, large uncertainties in 
the gas spin temperature, as well as non-uniformities 
within a single interstellar complex. The issue of the 
A^(II i) determination will be further discussed in Sec- 
tion [5] 

In Fig. [T] we show the differential emissivity per 
Hi atom measured in the Gould Belt. The results 
are compared with the predictions by GALPROP, a 
widely used code for CR propagation 0, Q . The re- 
sults are consistent with an independent study of LAT 
measurements at intermediate Galactic latitudes in 
the third quadrant Q • The H i emissivity in the Gould 
Belt obviously exceeds the GALPROP prediction by 
- 50%. 

However, the normalization of the spectrum is con- 
sistent with the expectations within the actual con- 
straints. In the energy range considered (200 MeV - 
10 GeV) the emission is dominated by the pionic com- 
ponent due to hadronic interactions. An important 
source of uncertainty comes from systematic errors in 
CR measurements: the current literature reports dif- 
ferences up to ^ 20% for the proton spectrum. An 
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Figure 1: Differential emissivity per Hi atom measured 
in the Gould Belt. Black points are our measurements, 
the shaded area gives the systematic uncertainties due to 
the LAT event selection efficiency. The lines represent 
the predictions by the GALPROP model 54_71Xvarh7S. 



additional ^ 30% uncertainty is due to the contribu- 
tion from heavier nuclei in both CRs and the ISM. The 
GALPROP model we adopted, following the method 
by Dermer d, @ , predicts an effective enhancement of 
'-^ 1.45 with respect to pure p-p emissivities. Different 
authors, however, report values as high as 1.75 — 2 [7]^. 

Provided a 50% scaling of the emissivity, the spec- 
tral shape is in good agreement with the expecta- 
tions. Thus we can conclude that CR nuclei in the 
local ISM are similar to those measured at Earth: 
further constraints will be given by improved direct 
measurements of CRs (by missions like PAMELA and 
AMS02). 

This is an additional evidence that LAT measure- 
ments are not consistent with the GeV excess seen by 
EGRET [see e.g. 1,0, as aheady deduced by LAT ob- 
servations at intermediate Galactic latitudes [l^l- A 
deviation with respect to the expected spectral shape 
of diffuse 7-ray emission was observed by EGRET over 
the whole sky for energies ^ 1 GeV, leading to sev- 
eral possible interpretations: instrumental effects, dif- 
ferences between typical CR spectra in the ISM and 
those in the proximities of the Earth, dark matter an- 
nihilation. The GeV excess was noticed also in the 
H I emissivity of local interstellar complexes [ll| , be- 
ing inconsistent with the results of our analysis. The 
same kind of difference is noticed in the spectra of 
bright point sources, notably for the Vela pulsar (T2j . 
and so the differences between the spectra measured 
by EGRET and those measured by the LAT are plau- 
sibly instrumental in origin. 



^The origin of the discrepancies is still unclear, either due 
to the calculation method (e.g. parametrization of the reaction 
yields versus numerical simulations) and/or to differences in the 
CR spectra within the experimental constraints. 



3. THE Xco GRADIENT IN THE OUTER 
GALAXY 

The molecular phase of the ISM is difffcult to trace, 
because its major constituent, H2, does not have ob- 
servable emission lines in its cold phase. The 2.6 mm 
line of the CO J = 1 — J> transition has often been 
used as a surrogate tracer of molecular masses. Obser- 
vations, notably virial mass calculations, suggest that 
the velocity- integrated antenna temperature Wqo j de- 
spite the optical thickness, is proportional to the total 
mass in the emitting region |13|. 

The conversion factor which transforms Wco inten- 
sity into 7V(H2) column density is known as Xqq = 
A^(H2)/VFco- For many years it has been considered 
as uniform across the Galaxy, but we have now ev- 
idence from virial masses [l4|, from COBE/DIRBE 
observations [l^, and from the Xco dependence on 
metallicity in external galaxies that it increases 
with Galactocentric radius. Such an increase was 
shown to be consistent with EGRET data [l3|, but 
the limited performance of previous 7-ray telescopes 
did not allow to get accurate measurements in the 
outer Galaxy [ll|, [3l ■ 

High-energy 7-ray observations provide important 
constraints on Xqq values. Since the molecular bind- 
ing energy is negligible with respect to the energy scale 
of the 7 radiation processes, the emissivity per H2 
molecule is twice the emissivity per H i atom for gas 
threaded by the same CR flux. This statement relies 
on the assumption that CR diffusion lengths are large 
enough to ensure similar CR spectra in the atomic and 
molecular phases of a complex (i.e. they do not vary 
significantly over distances of tens of pc) and that CRs 
penetrate molecular clouds uniformly to their core in 
spite of the enhanced magnetic fields (idea still under 
debate) . 

The results of our work allowed us to correlate the 
integrated emissivity per Wco unit {qco) and per Hi 
atom ((?Hi) in different regions along the line of sight, 
for several energy bands from 200 MeV to 10 GeV. 
The results are shown for the local arm in Fig. [5] 

The good linear correlation supports the assump- 
tion of similar spectra in the atomic and molecular 
phases of the clouds and the energy-independent pen- 
etration of CRs to the molecular cores. The slope of 
the linear relation provides an estimate of the Xqq 
ratio. Good linear correlations are also found for the 
Gould Belt clouds and in the Perseus arm. The val- 
ues of Xco in the different regions are summarized in 
Fig. [3] as a ftmction of Galactocentric radius. 

The results confirm a significant but moderate in- 
crease of Xco with Galactocentric radius. The gradi- 
ent is consistent with the metallicity decrease in the 
outer Milky Way and the strong relation seen between 
metallicity and Xco in external galaxies |l^. This 
relation is plausibly due to poor self-shielding against 
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Figure 2: Correlation between qco and qhi over 5 energy 
ranges between 200 MeV and 10 GeV in the local arm 
clouds. The solid red line represents the best linear fit. 
The fit uses the measurement uncertainties on both axes. 
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Figure 3: Black points: Xco values measured in the 
Gould belt, local and Perseus arm (from left to right). 
Blue dashed line: model by Strong et al. [13]. Green solid 
line: the conversion function based on virial masses by 
Nakanishi and Sofue UM- 



CO photodissociation by UV radiation in low metal- 
licity environments. 

Our results do not support the extremely large gra- 
dient proposed by Strong et al. beyond the solar 
circle (dashed blue line in Fig. ^ . The latter was in- 
troduced on the basis of non-7-ray data to alleviate 
the CR gradient problem (see below Section 2]). We 
find a value significantly lower than the prediction of 
this model for the molecular complex associated with 
NGC 7538 and Gas A in the Perseus arm, which is 
the most massive molecular complex beyond the solar 
circle. 

The values of Xco we measured are systematically 
lower than the conversion law derived by Nakanishi 
and Sofue [HI from virial masses. We have verified 
that a difference of comparable magnitude is found for 
our clouds between the masses derived from the Xqo 
ratios just measured and the virial masses (based on 
the same CO data used for the 7-ray analysis) [![. 



An instrumental origin of the discrepancy (e.g. cali- 
bration differences between different CO datasets) is 
therefore unlikely. As already said, the 7-ray masses 
may be biased by the non-uniform penetration of CRs 
within the dense cores of the clouds. The limited spa- 
tial resolution of the 7-ray data may also lead to sep- 
aration problems between the molecular phase of a 
cloud and its dense atomic envelope. On the other 
hand, virial masses rely on rather crude hypotheses: 
spherical clouds with simple density profiles and tur- 
bulent support to hold a cloud in equilibrium against 
gravitational collapse. The effects of magnetic sup- 
port are rarely considered (they were not in our cal- 
culation) and intrinsic velocity gradients can broaden 
the velocity dispersion. The virial estimate also cor- 
responds to the total dynamical mass, whereas the 
7-ray estimate corresponds to the mass spatially as- 
sociated with CO. Gas not properly accounted for by 
Wco intensities could explain, at least in part, the 
differences between the virial and 7-ray masses (see 
below Section O. 



4. THE COSMIC-RAY DENSITY 
GRADIENT IN THE OUTER GALAXY 

The CR origin is still unclear. The EGRET non- 
detection of the Small Magellanic Cloud [l^l, an ex- 
ternal galaxy member of the local group, highlighted 
that its CR content significantly differs from that of 
the Milky Way. The EGRET upper limit was recog- 
nized as evidence that CRs, at least below ~ 10^^ eV, 
have a Galactic origin. Supernova remnants (SNRs) 
have been considered for many years the best candi- 
dates as CR accelerators in our Galaxy. Whereas dif- 
fusive acceleration of electrons by SNR shock waves is 
obvious from multiwavelength observations, only re- 
cently we started accumulating indirect evidence for 
the acceleration of CR nuclei 7-ray observations in the 
TeV [HHP and GeV domain (24[. 

On the other hand, the distribution of SNRs in the 
Galaxy is very poorly determined because of the 
sparsity of the known sample (a few tens) and its large 
selection effects. Since the COS-B era it is well known 
that the gas 7-ray emissivity gradient across the outer 
Galaxy is much flatter than expected from the decline 
in SNR detections j26l | . 

Strong et al. [T3| proposed to use the observed 
pulsar distribution to trace CR sources because of 
their richer sample (27j . Large uncertainties in the 
distance derivation from the dispersion measurements 
can, however, easily bias their radial distribution. The 
decline in pulsar counts in the outer Galaxy is even 
steeper than the SNR one, strengthening the CR gra- 
dient problem. Strong et al. f7[ attempted to alleviate 
this problem by increasing the amount of molecular 
gas in the outer Galaxy, invoking a strong increase of 
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Xqo beyond the solar circle. Yet, tuning the Xqo 
ratio does not explain the flat qn i emissivity gradient 
recorded in the ubiquitous and massive atomic phase 
[26l nil ■ As discussed in Section [3] this solution is 
disfavored by LAT measurements. 

Our measurements of the Hi emissivity across the 
outer Galaxy are shown in Fig. 21 
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Figure 4: H i emissivities integrated between 200 MeV 
and 10 GeV (gui) measured in the Gould belt, local, 
Perseus and outer arm (from left to right). Black points 
are our measurements, the grey shaded area represents 
the systematic uncertainties in the LAT event selection 
efficiency. The blue dashed line gives the prediction by 
the GALPROP model 54_71Xvarh7S. 



The results are compared with the predictions by a 
GALPROP model which uses a CR source distribu- 
tion derived from pulsars. The CR gradient problem 
is confirmed by LAT data: the measured decline in 
emissivity is much flatter than the prediction by GAL- 
PROP for the pulsar (Fig. 15) and also the SNR dis- 
tribution (Fermi LAT collaboration, in preparation). 

Systematically underestimating the remote A^(Hi) 
densities as the result of self- absorption phenomena 
can lead to a bias in the emissivity gradient. Correct- 
ing for H I self absorption, Gibson et al. [13 suggested 
that a significant amount of H i mass has been over- 
looked by emission surveys in the second quadrant for 
the outer Galaxy. The large emissivity in the outer 
Galaxy may be explained also by a population of un- 
resolved sources clustering in the Perseus arm struc- 
tures. 

On the other hand, if the flat gradient is real, it 
challenges either the SNRs as GR parent population 
or the modeling of GR diffusion across the Galaxy. A 
possible explanation is provided by the large uncer- 
tainties on the SNR distribution. Additionally, the 
GR propagation parameters are usually derived from 
the isotopic composition of GRs measured at Earth 
[see e.g.|2|], a reasonable assumption, but not required 
by observational constraints. An alternative scenario 
is given by non-uniform CR diffusion fsol] . 



5. THE TRACERS OF THE 
INTERSTELLAR MEDIUM IN THE GOULD 
BELT 

The modeling of the 7-ray interstellar emission 
based on the combination of Hi and CO data is being 
challenged by observations. EGRET data, compared 
with other tracers, already suggested that a signifi- 
cant amount of gas is not properly traced by these 
radio/microwave lines [3l| . The nearby, off-plane, 
clouds located in the Gould Belt are well suited to 
probe ISM tracers, both because there is little con- 
fusion along the line of sight and because the higher 
linear resolution of the 7-ray data within the clouds 
allows a good separation of the emission arising from 
the different phases and sources. 

Following the method proposed by Grenier et 
al. [3l| , we have used interstellar dust as an additional 
gas tracer under the assumption that dust grains are 
well mixed with gas in the cold and warm phases of 
the ISM under study. We have used the total dust col- 
umn densities provided by the color excess E(B — V) 
map by Schlegel et al. [32]. From this map we have 
subtracted the parts linearly correlated with the best- 
fit combination of iV(Hi) and Wqo maps built for the 
four regions separated along the lines of sight, from 
the Gould Belt to the outer arm ^ij. 

The residual map is shown in Fig. [5] on the left. 
Within a few degrees from the Galactic plane, the 
map exhibits clumps of positive and negative resid- 
uals. Their origin is unclear because of the intrinsic 
limitations of the E(B — V) map near the plane [32| . 
At intermediate latitudes {\b\ > 5°), the map is dom- 
inated by structured envelopes of positive residuals 
surrounding the Gould Belt CO clouds. For compar- 
ison the Wco rnap is shown on the right of Fig. [5] 
Smaller negative residuals are found toward the dense 
CO cores. 

Adding the E(B — V) residual map significantly im- 
proves the fit to the 7-ray data Given the spatial 
correlation, we compared the 7-ray emissivities per 
unit of E(B — V) residuals (qebv) with the emissivi- 
ties per Hi atom (qui) in the Gould-Belt clouds. The 
results are shown in Fig. [51 

The good linear correlation between gEBV and ^hi 
over 2 decades in energy suggests that the same pro- 
cess produces the 7 rays associated with Hi and the 
dust-traced component. 

The 7-ray emission arising from hadronic interac- 
tions of CRs with dust grains is too faint to be de- 
tected. The IG emission coming from dust infrared 
emission up scattered by CR electrons is also expected 
to be faint and spectrally different from the 7r°-decay 
radiation associated with Hi 'sT]. Therefore, the re- 
sults point to the presence of additional gas which is 
not properly accounted for in the A^(Hi) and Wco 
maps. The negative residuals seen towards the CO 
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Figure 5: Maps of the region of Cassiopeia and Cepheus. Left: E(B — V) residuals (magnitudes), obtained from the 
map by Schlegel et al. [32| after the subtraction of the parts correlated with -N"(H i) and Wco maps HI. Right: Wco (K 
km s"'^) in the Gould Belt region, from the 2.6 mm survey of the CfA telescope [3^ . 
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Figure 6: Correlation between qebv and qm in the 
Gould Belt over 5 energy ranges between 200 MeV and 
10 GeV. The solid red line represents the best linear fit. 



cores can be interpreted as small local variations in 
the dust-to-gas ratio or in the dust thermal emissiv- 
ity in the denser, well-shielded parts of the molecular 
clouds. The nature of the additional gas, however, is 
not determined yet. 

Substantial uncertainties result from the approxi- 
mations applied to handle the radiative transport of 
the H I lines. The choice of a uniform spin temper- 
ature, Ts — 125 K, has often been adopted in the 
past, as in our study. Recent absorption surveys (sj] 
propose higher values in the outer Galaxy (250 — 400 
K) which would yield iV(Hi) densities lower than the 
present estimates. Conversely, colder and denser H i 
gas with Ts as low as 30 — 40 K is quite abundant in 



the mid- latitude clouds [35|. This phase contributes 
much larger TV (Hi) column-densities than the esti- 
mates based on a uniform spin temperature of a few 
hundreds K. Self absorption could also bias the A^(H i) 
column densities to low values. 

On the other hand, at the interface between the 
atomic and CO-traced phase of an interstellar com- 
plex we might find an additional phase where H2 ex- 
ists but it is not well mixed with CO (e.g. because 
H2 is more efficiently self-shielded than CO in the low 
extinction range) or densities are too low to excite 
the CO 2.6 mm line transition [s^. A poor correla- 
tion with Wco maps and, conversely, a good correla- 
tion with the E(B — V) color excess has already been 
shown for alternative molecular tracers (like the 9 cm 
line of CH) in translucent molecular clouds, and pos- 
sibly proposed for the translucent envelopes of giant 
molecular clouds [37| . 
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